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Abstract: We investigate the relation between a postulated skeleton expansion and the 
ji^i conformal limit of QCD. We begin by developing some consequences of an Abelian-like 

D I skeleton expansion, which allows one to disentangle running-coupling effects from the re- 

' maining skeleton coefficients. The latter are by construction renormalon-free, and hence 

hopefully better behaved. We consider a simple ansatz for the expansion, where an ob- 
^ ■ servable is written as a sum of integrals over the running-coupling. We show that in this 

framework one can set a unique Brodsky-Lepage-Mackenzie (ELM) scale-setting procedure 
as an approximation to the running-coupling integrals, where the BLM coefficients coincide 
with the skeleton ones. Alternatively, the running-coupling integrals can be approximated 
using the effective charge method. We discuss the limitations in disentangling running 
coupling effects in the absence of a diagrammatic construction of the skeleton expansion. 
Independently of the assumed skeleton structure we show that BLM coefficients coincide 
with the conformal coefficients defined in the small Pq (Banks-Zaks) limit where a pertur- 
bative infrared fixed-point is present. This interpretation of the BLM coefficients should 
explain their previously observed simplicity and smallness. Numerical examples are criti- 
cally discussed. 
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1 Introduction 



The large-order behavior of a perturbative expansion in gauge theories is inevitably 
dominated by the factorial growth of renormalon diagrams |jl], ^ 0, ^. In the case 
of quantum chromodynamics (QCD), the coefficients of perturbative expansions in 
the QCD coupling as can increase dramatically even at low orders. This fact, to- 
gether with the apparent freedom in the choice of renormalization scheme and renor- 
malization scales, limits the predictive power of perturbative calculations, even in 
applications involving large momentum transfer where as is effectively small. 

A number of theoretical approaches have been developed to reorganize the per- 
turbative expansions in an effort to improve the predictability of perturbative QCD. 
For example, optimized scale and scheme choices have been proposed, such as the 
method of effective charges [ECH] [Q, the principle of minimal sensitivity [PMS] 0, 
and the Brodsky-Lepage-Mackenzie [BLM] scale-setting prescription and its gen- 
eralizations P|-pO[. More recent developments |^ include the resummation of the 
formally divergent renormalon series and the parameterization of related higher-twist 
power-suppressed contributions. 

In general, a factorially divergent renormalon series arises when one integrates over 
the logarithmically running coupling as{k'^) in a loop diagram. Such contributions 
do not occur in conformally invariant theories which have a constant coupling. Of 
course, in the physical theory, the QCD coupling does run. Nevertheless, relying on 
a postulated "dressed skeleton expansion", we shall show that a conformal series is 
directly relevant to physical QCD predictions. 

In quantum electrodynamics the dressed skeleton expansion can replace the stan- 
dard perturbative expansion. The skeleton diagrams are defined as those Feynman 
graphs where the three-point vertex and the lepton and photon propagators have 



no substructure |21|]. Thanks to the QED Ward identity, the renormalization of the 
vertex cancels against the lepton self-energy, while the effect of dressing the photons 
in the skeleton diagrams by vacuum polarization insertions can be computed by inte- 
grating over the Gell-Mann Low effective charge a{k'^). The perturbative coefficients 
defined from the skeleton graphs themselves are conformal - they correspond to the 
series in a theory with a zero (3 function. Therefore they are entirely free of running 
coupling effects such as renormalons. Each term in the dressed skeleton expansion 
resums renormalon diagrams to all orders in a renormalization scheme invariant way. 
The resummation ambiguity, which is associated with scales where the coupling be- 
comes strong, can be resolved only at the non-perturbative level. 

In QCD, a skeleton expansion can presumably be constructed based on several 



different dressed Green functions (see [^). A much more interesting possibility, 



which is yet speculative, is the existence of an Abelian-like skeleton expansion, with 
only one effective charge function. The construction of such an expansion is not 
straightforward due to the presence of gluon self-interaction diagrams and the essen- 
tial difference between vacuum polarization insertions and charge renormalization. 
Nevertheless, at the one-loop level there is a diagrammatic algorithm, the so-called 
"pinch technique" [p3| , that allows one to identify in every non-Abelian diagram the 
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part which can be absorbed into the renormahzation of the effective gluon propagator. 
The sum of all the vacuum-polarization-like parts turns out to be gauge-invariant, 
thus defining a natural candidate for the non-Abelian equivalent of the Gell-Mann 
Low effective charge, q:s(/c^). Moreover, the pinch technique leads to Ward identities 
similar to those of the Abelian theory: after the vacuum-polarization-like parts have 
been taken into account, the remaining vertex correction cancels against the quark 
self-energy. In this way the pinch technique achieves the first step in the construction 
of an Abelian-like skeleton expansion. Recently there have been some encouraging 
developments [0, in the application of the pinch technique beyond one-loop and 
its possible relation to the background field method. The hope is that these tech- 
niques will eventually provide a proof of existence for the skeleton expansion as well 
as an all-order constructive definition for the non-Abelian skeleton structure and the 
non-Abelian skeleton effective charge as(fc^). 

In this paper, we shall postulate that an Abelian-like skeleton expansion can be 
defined at arbitrary order in QCD. We shall not deal here with the diagrammatic 
construction of the skeleton expansion but rather restrict ourselves to the conse- 
quences which follow from such a structure. To this end we will introduce a simple 
ansatz for the skeleton expansion, where similarly to the Abelian case, a generic ob- 
servable is written as a sum of integrals over the running coupling. As in QED, we 
can then identify running coupling effects to all orders, and treat them separately 
from the conformal part of the perturbative expansion. A considerable simplification 
is achieved, for instance, by assuming that the dependence on the number of light 
quark flavors Nf originates only in the running coupling itself, as in the Abelian the- 
ory with light-by- light diagrams being excluded. As a consequence, the coefficients 
appearing in the assumed skeleton expansion are Nf independent. By construction 
these skeleton coefficients are free of renormalons and are therefore expected to be 
better behaved. We will show that they have a simple interpretation in the presence 
of a perturbative infrared fixed-point, as occurs in the small (3q limit: they are the 
"conformal" coefficients in the series relating the fixed-point value of the observable 
under consideration with that of the skeleton effective charge. Thus, given the as- 
sumption that these coefficients are Nf independent, they can be obtained from the 
standard perturbative coefficients using the Banks-Zaks expansion [^, ^ where 
the fixed-point coupling is expanded in powers of [3q. 



The conformal series can be seen as a template |T0[ for physical QCD predic- 
tions, where instead of the fixed coupling one has at each order a weighted average 
of the skeleton effective charge as{k'^) with respect to an observable- (and order-) 
dependent momentum distribution function. The momentum integral corresponding 
to each skeleton term is renormalization-scheme invariant. Had the skeleton effec- 
tive charge been known at all scales, this integral could have been unambiguously 
evaluated, thus including both perturbative and non-perturbative contributions. In 
practice it can be evaluated up to power-suppressed ambiguities, which originate in 
the infrared where the coupling becomes strong. These infrared renormalon am- 
biguities can be resolved only by explicitly taking non-perturbative effects into ac- 
count. Since such effects cannot be calculated with present methods, they can only 
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be parametrized. Indeed, a natural parametrization in the form of an infrared finite 
coupling emerges from the structure of the skeleton integral. This way the skele- 
ton expansion gives a natural framework in which renormalon resummation and the 



analysis of non-perturbative power corrections are performed together |31]. 

As an alternative to computing a dressed skeleton integral, one can approximate 
it by evaluating the coupling at the BLM scale , in analogy to the mean- value theo- 
rem [11]. By going to higher orders in the perturbative expansion, this approximation 
can be systematically improved. Another possibility considered here is to approxi- 
mate each dressed skeleton term separately using the effective charge approach. This 
approach is tailored to deal with running coupling effects and it by-passes remain- 
ing scheme and scale setting ambiguities in the power series expressions for the BLM 
scales. Assuming a simple form of the skeleton expansion, running coupling effects 
can be disentangled from the remaining conformal expansion by tracing the Nj de- 
pendence of the coefficients. In this case BLM scale-setting (or the ECH alternative) 
can be applied to a generic QCD observable based on the knowledge of the first few 
coefficients. In the general case disentangling running coupling effects becomes more 
involved, and it eventually requires a diagrammatic construction of the skeleton ex- 
pansion. We emphasize that both the BLM scale-setting method and the suggested 
ECH method remain on the perturbative level and, as opposed to the infrared finite 
coupling approach mentioned above, these methods are not particularly suited to 
deal with renormalon ambiguities and the related power-corrections. 

BLM scale-setting can also be applied to the perturbative relation between the 
effective charges of two physical observables. This results in a specific "commensu- 
rate scale relation" between the two quantities. The coefficients appearing in 
such relations are conformal and, as guaranteed by the transitivity property of the 
renormalization group, they do not depend on the intermediate scheme used. This 
way conformal relations appear to be relevant for real-world QCD predictions even 
in the absence of complete understanding of the underlying skeleton structure. In 
the case of the Crewther relation [P2| , |TB|, which connects the effective charges of 
the e~^e~ annihilation cross section to the Bjorken and Gross-Llewellyn Smith sum 
rules for deep inelastic scattering, the conformal relation is simply a geometric series. 
This example highlights the power of characterizing QCD perturbative expansions in 
terms of conformal coefficients. 

This paper is organized as follows: we begin in section 2 by recalling the concept 
of the skeleton expansion in the Abelian case |2T| and stating the main assumptions 
concerning the non- Abelian case. We continue, in section 3, by reviewing the standard 
BLM scale-setting procedure and recalling the ambiguity of the procedure beyond the 
next-to-leading order. We then show how this ambiguity is resolved upon assuming a 
skeleton expansion, provided we work in the appropriate renormalization scheme, the 
"skeleton scheme" , and require a one-to-one correspondence between the terms in the 
BLM series and the dressed skeletons. We also discuss in this section the limitations 
of applying the formal BLM procedure in the absence of a diagrammatic construction 
of the skeleton expansion. 

In section 4 we present an alternative to performing an explicit scale setting. 
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by using the ECH method as a tool to resum running-couphng effects within the 
framework of the assumed skeleton expansion. Close connections between the two 
approaches are pointed out. In Appendix A we look at the original ECH approach 
from the point of view of the skeleton expansion, comparing it to the application of 
the ECH method to the leading skeleton term. We also calculate there the /3o = 
limit of the skeleton coupling (3 function coefficient (32- 

We then come to the main subject of the paper, the relation between the coef- 
ficients which remain after applying BLM scale-setting and the conformal limit of 
QCD. We derive (section 5) a relation between these BLM coefficients and the con- 
formal coefficients defined in the infrared limit in the conformal window, where a 
non-trivial perturbative fixed-point exists In section 6 we show explicitly 

that the conformal coefficients, calculated using the Banks-Zaks expansion, are the 
same as the ones in the BLM series. In section 7 we recall previous observations con- 
cerning the smallness of conformal and Banks-Zaks coefficients, and examine whether 
this apparent convergence can be explained by the absence of renormalons in such 
relations. The conclusions are given in section 8. 



2 Renormalons and the skeleton expansion 

Consider a Euclidean QED observable a^((5^), which depends on a single external 
space-like momentum and is normalized as an effective charge. The perturbative 
expansion in a generic renormalization scheme is then given by, 

= + ria{^^f + r2a{fi'^f H , (1) 

where a = a/n and fi is the renormalization scale. 

The perturbative series can be reorganized and written in the form of a skeleton 
expansion 

aniQ') = RoiQ') + sMQ') + S2R2iQ') + ■ ■ ■ , (2) 

where the first term, i^o^ corresponds to a single dressed photon: it is the infinite 
set of "renormalon diagrams" obtained by all possible vacuum polarization insertions 
into a single photon line. The second term, siRi, corresponds to a double dressed- 
photon exchange and so on. In QED, vacuum polarization insertions amount to 
charge renormalization. Thus Rq can be written as 



Rom^ f-^ik')4>o{kVQ')^ (3) 



where k"^ is the virtuality of the exchanged photon, a{k'^) is the Gell-Mann Low 
effective charge representing the full propagator, and 0o is the (observable dependent) 
Feynman integrand for a single photon exchange diagram, which is interpreted as the 
photon momentum distribution function Similarly, Ri is given by 



R,iQ') ^ f -aikf) -a{kl) {k\IQ\ kl/Q^) ^ ^ (4) 



4 



and so on. 

For convenience the normalization of 0j in Ri{Q^) has been set to 1 such that 
the Ri{Q^) in have an expansion Ri{Q'^) = a{Q'^y~^^ + ■ ■ ■• For example, the 
normalization of 0o(^^/Q^) in Ro is 

„ >o(^VQ^)^ = l. (5) 

In QED fermion loops appear either dressing the exchanged photons or in light- 
by-light type diagrams, where they are attached to four or more photons (an even 
number). Barring the latter, the dependence on the number of massless fermion 
flavors Nf is fully contained in the Gell-Mann Low effective charge. It follows that 
the skeleton coefficients Si as well as the momentum distribution functions (pi are 
entirely free of Nf dependence. Light-by-light type diagrams have to be treated 
separately, as the starting point of new skeleton structures. 

The skeleton expansion is a renormalization group invariant expansion: each 
term is by itself scheme invariant. This is in contrast with the standard scale and 
scheme dependent perturbative expansion (0). The renormalons in (|l]) can be ob- 
tained upon expansion of the dressed skeleton terms in in some scheme. Let us 
consider first the leading skeleton and examine, for simplicity, its expansion in 
a((5^). We assume that a(/c^) obeys the renormalization group equation. 



da{k 



d\nk^ -(/3oa(A:')' + Aa(A;=^)' + «A;2)4 + ...) = p(a) (6) 
where /3o is negative in QED and positive in QCD. Then a{k'^) can be expanded as 

a{k^) = a{Q^)+(3ota{Qy + {f3it + /^qV) a{Q^)''+ (^^st + ^A/So^' + /^o^') d{QY+- ' ' 

(7) 

where t = — ln(/c^/Q^). Inserting this in eq. (^) under the integration sign we obtain 
RoiQ') = d{Q') + ri'^Podm' + (rf + r^'^p,) a^Q^ 

+ {rlt^P'o + lr?PiPo + r^^^P^) + ■■ (8) 



where 



' i 



2 



\n(kyQ')r Mk'm^. (9) 



fc2 



We note that in the large Nf (large /3o) limit[|, the perturbative coefficients ■ 



Si) 



*In QCD, Abelian correspondence in the large Nf limit requires that the coefficient Pi of the 
skeleton coupling (3 function would not contain Nj'^^. It has to be a polynomial of order Nj 
in Nf. This would guarantee that in the large Nf limit /3{a) is just the one-loop f3 function. Note 
that while some schemes (e.g. MS and static potential effective charge) have this property, generic 
effective charges (defined through observable quantities) do not. This property of the skeleton 
scheme is used making the identification of r|''' in (0) as the large Nf coefficients. 
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and thu^ 



large /3o 



Li=0 



(10) 



At large orders i ^ 1, both small and large momentum regions become dominant 
in (P), giving rise to the characteristic renormalon factorial divergence {rf' ~ 
As mentioned above, this is believed to be the dominant source of divergence of the 
perturbative expansion (|l]). On the other hand, in the skeleton expansion (H) the 
renormalons are by definition resummed and so the remaining coefficients Si should 
be free of this divergence. These coefficients are expected to increase much slower 
leading to a better behaved expansion. 

As mentioned in the introduction, the generalization of the Abelian skeleton ex- 
pansion to QCD is not straightforward. Diagrammatically, the skeleton expansion in 
QCD has a simple realization only in the large Nf limit where gluon self-interaction 
contributions are negligible so that the theory resembles QED|^ In the framework 
of renormalon calculus, one returns from the large Nf limit to real world QCD by 
replacing Nf with the linear combination of Nf and Ca = N^ which appears in the 
leading coefficient EDI of the (3 function. 



/?o = ^(fc^-^iV/). (11) 



This replacement, usually called "naive non-Abelianization" |T^, 0, amounts 
to taking into account a gauge invariant set of diagrams which is responsible for the 
one- loop running of the coupling constant. 

To go beyond the "naive non-Abelianization" level constructing an Abelian-like 
skeleton expansion in QCD, one needs a method to identify skeleton structures and 
to isolate vacuum-polarization-like insertions which are responsible for the running of 
the coupling at any order. The pinch technique |2^, ^ ^ may provide a systematic 
way to make this identification. The resulting set of skeleton structures would surely 
be larger than in the Abelian theory. It may include, for example, fermion loops at- 
tached to an odd number of gluons, which vanish in the Abelian limit. Like Abelian 
light-by-light type diagrams, these structures should be treated separately. As op- 
posed to the Abelian theory, where light-by-light type diagrams are distinguished by 
their characteristic dependence on the charges, in the non- Abelian case these struc- 
tures may not be separable based only on their color group structure. We assume 
that there is a unique (gauge invariant) way to identify skeleton structures in QCD, 
making eq. (|^) relevant. In much of the discussion that follows we shall further make 
the assumption that the entire dependence of aji{Q^) on Nf is through the running 
coupling. Thus in our ansatz Sj and (pi are Nf independent, just like in the Abelian 
case with light-by-light diagrams being excluded. Of course, the class of diagrams 



'I' We comment that the sub-leading terms in in (^) of the form /3i/3q~^ were computed to 

all-orders in |p^ . However, other terms which involve higher order coefficients of the /3 function 
contribute at the same level in 1/Nf. 

tThis can also be understood from the Ca limit discussed in ref. [p9[. 
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containing fermion loops as part of the skeleton structure should eventually be taken 
into account. 

Another simplifying assumption we made already in the Abelian case is that our 
ansatz (0) contains only one skeleton at each order, whereas in general there will be 
several skeletons contributing at each order. The simplest example is e~e~ scattering 
with both t- and w-channel exchange. Several skeletons at the same order also occur 
in single-scale observables considered here, and therefore (^) should be generalized 
accordingly. We shall return to this point in the next section. 

We stress that the coupling constant a(/c^) in (|^) is understood to be a specific 
effective charge, in analogy to the Gell-Mann Low effective charge in QED. This 
"skeleton effective charge" a(/c^) should be defined diagrammatically order by order 
in perturbation theory. In the framework of the pinch technique, a(fc^) has been 
identified at the one-loop levelfl e.g. it is related to the MS coupling by 



/ e 5\ Ca 



Ams (/i')' + ■ ■ • (12) 



Recently, there have been encouraging developments ^ in the application of the 
pinch technique beyond one-loop. This would hopefully lead to a systematic identifi- 
cation of the "skeleton effective charge" at higher orders, namely the determination 
of higher order coefficients {f3i for i > 2) of the f3 function (3{a) = da/dlnk^. This f3 
function should coincide with the Gell-Mann Low function upon taking the Abelian 
limit Ca = (see ref. 

Being scheme invariant and free of renormalon divergence, the skeleton expansion 
(0) seems much favorable over the standard perturbative QCD expansion (|I]). This 
advantage may become crucial in certain applications, e.g. for the extraction of as 
from event shape variables |PT|. However, in the absence of a concrete all-order di- 



agrammatic definition for the skeleton expansion in QCD, running-coupling effects 
cannot be systematically resummed beyond the single dressed gluon level correspond- 
ing to the leading skeleton. In particular, the momentum distribution functions of 
the sub-leading skeletons are not known. On the other hand, the BLM scale-setting 
procedure, which is well defined at higher orders, can be considered as a manifesta- 
tion of the skeleton expansion. As we shall see, it is possible in this framework to 
approximate the sub-leading skeleton terms, provided the correct skeleton scheme is 
used. Currently, since the skeleton effective charge has not been identified, the choice 
of scheme in the BLM procedure remains an additional essential ingredient. 



3 BLM scale-setting 

The BLM approach is motivated by the skeleton expansion. The basic idea is 
that the dressed skeleton integral (^ can be well approximated by Rq — a(yU^) + ■ ■ ■ 
provided that the renormalization scale fi is properly chosen. Indeed, by the mean 

§This means that the corresponding QCD scale A is identified. 
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value theorem [|TT1, there exists a scale ko such that 



RoiQ') = f -a{e) <Po {kVQ') ^ = -aikl) f M^VQ')^ = a{kl) (13) 

where the last step follows from the assumed normalization for 0j (^. 

A first approximation to /cq is given by the average virtuality of the exchanged 
gluon, 

kl, = Q'e^V (/f l^^'^o(W)^ /^°°0o(W)^) =g^exp(-r«) 

(14) 

where is the next-to- leading coefficient of aR in the large (3q limit (P). The scale 
([T^ ) is called the "leading order BLM scale" . It can be determined directly from the 
Nf dependent part of the next-to- leading coefficient (ri) in the perturbative series of 
the observable in terms of a{Q'^), 

aR{Q') = a(Q2) + r^aiQ^f + r^a^Q^f + ■■■. (15) 

Thanks to the linear Nf dependence of ri, it can be uniquely decomposed into a term 
linear in (3q, which is related to the leading skeleton, and a free term 

ri = ri°) + ri^)/3o, (16) 

where both r^^ and rf'^ are Nf independent. After BLM scale-setting, with /cqq given 
by (0), one has 

aniQ') = -a{kl,) + r?h{kl,f + ■ ■ ■ . (17) 

Thus, technically, the BLM scale-setting procedure amounts, at leading order, to 
eliminating the /3o dependent part from the next-to-leading order coefficient. Note 
that although the leading order BLM scale /co,o of (li) h^^s a precise meaning as the 
average gluon virtuality, it is just the lowest order approximation to ko of eq. (|TB|). 
In other words, aiming at the evaluation of the leading skeleton term (^, setting the 
scale as fco,o is just the first step. This approximation can be systematically improved 
(see eq. (|33| ) below) in higher orders. 

3.1 Multi-scale BLM and skeleton expansion correspondence 



A BLM series |T^ can be written up to arbitrary high order 



aij(Q') = a{kl) + cia{kl) + C2a{kl) + c^a{kl) + ■ ■ ■ (18) 

where the kf are, in general, different scales proportional to the external scale (as in 
([T^)) and Cj are Nf independent coefficients. The intuition behind this generalization 
is that each skeleton term in (|^) is approximated by a corresponding term in the multi- 
scale BLM series: each skeleton term may have different characteristic momenta. This 
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one-to-one correspondence with the skeleton expansion requires that the couphng a 
will be the skeleton effective charge a = a such that 

R,{Q^) ^ -a{e, y+\ (19) 

In this case the coefficients of sub-leading terms in (|18D should coincide with the 
coefficients of the sub-leading skeleton terms, namely q = s^. 

More generally, a BLM series can be formally written in an arbitrary scheme: 
then the coupling a in (|l^) can be either defined in a standard scheme like MS or, 
as suggested in |T^, be another measurable effective charge. However, in such cases 
there is no direct correspondence with the skeleton expansion (Q) , and as a result the 
forthcoming motivation for a unique scale setting is lost. 

Let us recall how the BLM scale-setting procedure is performed beyond the next- 
to- leading order ||T^, ^, yielding an expansion of the form (^). Suppose that the 
perturbative expansion of a/j(Q^) in terms of a(Q^) is given byp] 

aR{Q^) = a{Q^) + na{Q^f + r2a{Q^f + r:,a{Q^f + ■ • • (20) 

Based on the fact that are polynomials of order i in Nj and that (3o and (3i are 
linear in iVj, we can write ri as in (|16D and 

r, = rf + r«rf )/?o + /^o' + r«A (21) 

where rp'' are Nf independent. The reason for the (3i dependent term in (^) shall 
become clear below. Expanding a{kf) in terms of a(Q^) similarly to eq. (|^), the 
next-to-next-to-leading order BLM series (p!8[) can be written as 

aniQ^) = a(Q2) + (^ + toPo) a(Q')' + (c2 + 2hciPo + toPi + ^P^) a(Q')l (22) 

Writing the scale-shifts ti = ln{Q'^/kf) as a power series in the coupling 

U = U^o + hi aiQ^) + ti,2 a(Q')' + ■ ■ ■ (23) 

where ti^ are assumed to be Nf independent, we get 

aniQ^) = a{Q^) + {ci + to,of3o)a{Qy (24) 
+ (c2 + (2ti,oCi + to,i)/5o + hoPi + to,o/3o) «(<5' 



-)2\3 



An order by order comparison of (|2^) and (^Dj) yields the scale shifts to = lii(Q^/^o) 
and ti = ln{Q'^/kl) and the coefficients Ci and C2 in terms of ri and r2 and the 
coefficients of the (3 function of a{Q'^). The comparison at the next-to-leading order 
gives 

ci = (25) 

^We work now in a generic scheme but in contrast to (|l|) we start here with the renormahzation 
scale fJ, = Q thereby simpHfying the formulas that follow. Since the scale is tuned in the BLM pro- 
cedure, this initial choice is of little significance. The only place where the arbitrary renormahzation 
scale is left at the end is in the power series for the scales-shifts, eq. ( p3|) below. 
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and 



tofl = r[". (26) 
The comparison at the next-to-next-to-leading order for the Pi independent piece 



gives 



C2 — / 2 

while for the /?o dependent piece it gives 



(0) 



(27) 



(28) 



Thanks to the explicit /3i dependent term introduced in (pi]), the equality of the 
corresponding piece there to that in ( p4D is satisfied based on the next-to-leading 
order result (|26|) . To proceed we need to specify to,i and ti^ such that eq. ( p8|) is 
satisfied. Having two free parameters with just one constraint there is apparently 
no unique solution. Two natural possibilities are the so called multi-scale BLM 



prescription [12 



^0,1 
^1,0 



1(1) 



— I r 



,(1) 



(29) 



and the single-scale BLM prescription M where ti^ = to,o and 



,(2) 



-2rfVf)+r«rf). 



(30) 



Having in mind the original motivation for BLM, it is interesting to examine the 
case where the scheme of a coincides with the skeleton effective charge d. Then 
we would like to have a one-to-one correspondence (|l^) between the terms in the 
BLM series (|18|) and those of the skeleton expansion (^. The multi-scale procedure 
is consistent with this requirement: the leading term a(fco) in the BLM series ([l^) 
represents only the leading skeleton term Rq in (EI), since the scale-shift 



+ 



,(2) 



(31) 



involves only coefficients which are leading in the large (3o limit and originate in 0o 
(cf. eq. (D). On the other hand the single-scale procedure violates this requirement, 
since there to involves (|30|) terms which are sub-leading in /5o and do not belong to 
the leading skeleton term i?o- In fact, in order to guarantee that the scale-shift to 
would represent just the leading skeleton Rq we are bound to choose to,i proportional 
to Po and thus the solution ( ^9]) is uniquely determined. 

We see that a unique scale-setting procedure at the next-to-next-to-leading order 
(r2) is implied by the requirement that the scale-shift to should represent the leading 
skeleton Ro- In order to continue and apply BLM at the next order [r^) we have to 
impose further constraints based on the structure of both Ro and -Ri. 
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3.2 BLM scale-setting for the leading skeleton 



Let us first examine the structure of tfie scale-sliift to by applying BLM to a liypo- 
tlietical observable tfiat contains only an Rq term of the form (^). Expanding the 
coupling a(A;^) under the integration sign in terms of a((5^) we obtain We would 
like to apply BLM to the latter series obtaining simply a(fco); with to = ln(Q^/fcQ) = 



^0,0 + ^o,i'^(*5^) + ■ ■ ■• Expanding a{kl) we obtain from 
D with (|3|) we get 



a{kl) = 
+ 

Comparing 

to = 
+ 



(32) 



,(3) 



^(2) 
'2 



(33) 



(i)V 



^(2) 
'2 



,»)y 



.ft a(Ql' + 



Here we recovered the two leading orders in to of eq. (0). At order a{Q^Y 

obtained an explicit dependence on both /5o and The combination rg^-* — [r'l''^ 
appearing at the next-to-leading order in to has an interpretation as the width of 
the distribution 0o, assuming the latter is positive definite (see [1^, |l3)- In general, 
eq. (|33| ) can be written in terms of central moments of the distribution 0O; defined by 




1 ISj , 



In 



00 / 



for n > 2, where Mi 
the central moments we have 



In^ 



00 



In p- corresponds to rj^'* in 



eq. 



(34) 



In terms of 



tn 



Ml + M2(3o aiQn + UMs + M^M^] (3^ + -M^Pi \ a(Q 



+ ■ ■ 



2\2 



+ ■ 



(35) 



where in the second step we changed the scale from to the leading order BLM 
scale /cg simpler expressions for the coefficients of the to series. At large orders 

n the moments M„ become sensitive to extremely large and small momenta and 
thus develop renormalon factorial divergence, similarly to the standard perturbative 
coefficients in eq. (|^). We thus see that in the BLM approach, the scale-shift itself is 
an asymptotic expansion, affected by renormalons. 



3.3 BLM scale-setting for sub-leading skeletons 

Next, let us consider an Ri term, given by (^). Expanding the couplings a(kl) and 
(3(^2) under the integral in terms of a(Q^) using (0), we get (cf. the expansion of Rq 
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in eq. (|)) 



where 



' 2 

^(2) 



^(1.0) ^ ^(0,1) 



(36) 



(37) 



with 



Hkl/Q')]' [~Hkl/Q')f 01 {kl/Q',kl/Q 



U2 

t^l 1^2 



(3^ 



The BLM scale-setting procedure can now be apphed according to (p^Of): RiiQ"^) 
given in eq. ( p6| ) should be written as aiklY . Expanding a{k\Y in terms of aiQ"^) 
using (1^) and ti = ti Q + ti,ia(Q^) + ■ ■ ■ we have 



-a{kiY = am' + 2ti,o/3oa(Q^)^ + 2ti,i/?o + 3tf o/?o + 2ti,oA «(Q ) + 



The comparison with (|36D at the next-to-leading order implies 



(1) 



^1,0 = 

The comparison at the next-to-next-to-leading order then yields 



(39) 



(40) 



(41) 



which implies that ti^i, just as to,i; is bound to be proportional to I3q. Finally we 
obtain the scale-shift for 



= -A> + - 



- \ ( 



(1) 



/3oa(Q'). 



Similarly, applying BLM to i?2, 

i?2 = a(Q')' + r?^/5oa(Q')' + 

we get 



1 



to = -To 

3 ^ 



(1) 



(42) 

(43) 
(44) 
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3.4 Skeleton decomposition and its limitations 



Let us now return to the case of a generic observable (|20|) and see that with these 
skeleton-expansion-correspondence constraints there is a unique BLM scale-setting 
procedure. The basic idea is that, given the existence of a skeleton expansion, it is 
possible to separate the entire series into terms which originate in specific skeleton 
terms. This corresponds to a specific decomposition of each perturbative coefficient 
Tj similarly to (0) and (0). Then the application of BLM to the separate skeleton 
terms, namely representing Ri by d{kfy~^^, immediately implies a specific BLM scale- 
setting procedure for the observable. For example, when this procedure is applied 
up to order a(Q^)^, the scale-shifts for i = 0,1,2 are given by (0), (|4^), 
respectively. 

To demonstrate this argument let us simply add up the expanded form of the 
skeleton terms up to order a{Q^)'^ with Rq given by (||), Ri by (^) and R2 by (|43|) . 
For i?3 we simply have at this order R3 = a(Q^)^. Altogether we obtain. 



a + 



+ 

+ 



.(1), 



a' 



52 + s,r'2>Po + r!^'>P'o+rr'Pi^ 

53 + S2r?Po + sirf /3o' + rfP! + r?P2 + PiPo + sA'^ Pi 



(45) 



.(1), 



Here we identify the notation Si which is the coefficient in front of the skeleton term 
Ri with rf \ We recognize the form of ri and r2 as the decompositions introduced 
before in eq. (|T6|) and pH ) in order to facilitate the application of BLM. We see 
that the skeleton expansion structure implies a specific decomposition. Suppose for 
example we know ri through in the skeleton scheme. Eq. (^) then defines a unique 
way to decompose them so that each term corresponds specifically to a given term in 
the skeleton expansion. The decomposition of rj includes a polynomial in /3o up to 
order /3q, 



fc=i 



where Sq = 1 by the assumed normalization. The other terms in 
explicitly on higher coefficients of the /3 function (3j with 1 < j < i - 



(46) 

in (^Sp depend 
1. Up to order 

a{Q'^)'^ these terms depend exclusivelyjU] on coefficients rj^'' which appeared at previous 
orders in the jSo polynomials (0). Finally, we need to verify that a decomposition of 
the form (|45|) is indeed possible. For a generic observable qr, the coefficient is a 
polynomial of order i in Nf. Since the j3 function coefficients Pi are also polynomials 
of maximal order i, the decomposition of according to ( ^5]) amounts to solving i + 1 
equations with i + 1 unknowns: r^'^^ 
solution. 



with < k < i. Thus in general there is a unique 



"As we shall see below, this is no longer true beyond this order, where the coefficients depend on 
moments which appeared at previous orders, but cannot be expressed in terms of the lower order 
coefficients themselves. 
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We see that based on the assumed skeleton structure, one can uniquely perform 
a "skeleton decomposition" and thus also BLM scale-setting which satisfies a one- 
to-one correspondence of the form ([T9|) with the skeleton terms. By construction in 
this procedure the scale to is determined exclusively by the large (3q terms which 
belong to Ro (see (p3D), ti is determined by rf"^"* terms which belong to Ri (see 
(|4^)), t2 is determined by rf"^-* terms which belong to R2, etc. 

It should be stressed that formally the decomposition (^S]), and thus also BLM 
scale-setting, can be performed in any scheme: given the coefficients up to order 
n, all the coefficients Sj and rj-''' for i < n and j < i are uniquely determined. No 
special properties of the "skeleton effective charge" were necessary to show that the 
decomposition is possible. Even the assumption that for this effective charge the 
f3 function coefficients f3i are polynomials of order i can be relaxed. For example, 
the decomposition (^51) can be formally performed in physical schemes where /3j are 



polynomials of order i + In this case, however, the interpretation of r|"'^ in terms of 
the log-moments of distribution functions is not straightforward. It is also clear that 
a one-to-one correspondence between BLM and the skeleton expansion (|19D exists 
only if the coupling a is chosen as the skeleton effective charge d. 

Let us now address several complications that limit the applicability of the above 
discussion. First, we recall the assumption we made that the entire dependence of the 
perturbative coefficients on Nj is related to the running coupling. This means that 
any explicit Nj dependence which is part of the skeleton structure is excluded from 
(|45|). In reality there may be skeletons with fermion loops as part of the structure, 
which would have to be identified and treated separately. 

Having excluded such Nf dependence, we have seen that up to order d{Q^y a for- 
mal "skeleton decomposition" (^51) of the perturbative coefficients can be performed 
algebraically without further diagrammatic identification of the skeleton structure. 
This is no longer true at order a(Q^)^, where the "skeleton decomposition" requires 
the moments of the momentum distribution functions to be identified separately. Such 
an identification depends on a diagrammatic understanding of the skeleton structure. 



Looking at -Ri, the coefficient of a{Q^)^ in eq. ( ^6]) is 



/5o^ [<P?''^ + <P?''^ + 0^'^ + <P?''^] (47) 



+ P2 



^(1,0) ^ ^(0,1) 



Writing the a(Q^)^ term in (|i5|), one will find as before, that the terms which depend 
explicitly on higher coefficients of the f3 function jSi with 1 < / < 3, contain only 
moments of the skeleton momentum distribution functions 0^'^^ which appeared in 
the decomposition (^) in the coefficients of /Sq'*''^ at the previous orders. 

However, the coefficient of /3i/3o will depend on a new linear combination of moments, 
different from the one identified at order a(Q^)^ (compare the coefficient of /3i/?o in 
(|47| ) with r^^ in eq. (0)). Thus, strictly based on the algebraic decomposition of the 



coefficients at previous orders there is no way to determine the coefficient of PiPo at 
order a(Q^)^. Additional information, namely the values of 4'i''^\ (f)'i'^^ and 01°'^'' is 
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required. In the Abelian case, where the diagrammatic identification of the skeleton 
structure is transparent, it should be straightforward to calculate these moments 
separately. In the non-Abelian theory this not yet achievable. 

The need to identify the skeleton structure, as a preliminary stage to writing the 
decomposition of the coefficients (and thus also to BLM scale-setting) may actually 
arise at lower orders if several skeletons appear at the same order. As mentioned 
in the previous section, even in the Abelian case the assumed form of the skeleton 
expansion (|^) is oversimplified in this sense and should be generalized to include 
several different SiRi{Q'^) terms at any order i. In the non-Abelian case one should 
expect the set of skeleton diagrams to be larger. 

The simplest possibility to imagine is that the momentum distribution functions 
(pi are N,, independent. It is then natural to expect that at any given order there 
will be several skeletons, where each of them is characterized by its own color group 
structure. For example, let us assume that in the case of the QCD correction to the 
photon vacuum polarization sii?i(Q^) should be replaced by a sum of three skeleton 
terms, SiRiiQ"^) + Si^RYiQ"^) + s\^R^^{Q'^), where the first two terms correspond to 
double gluon exchange (which exist in the Abelian theory) - the planar (p) and the 
non-planar (np) skeleton diagrams, and the last term corresponds to the three gluon 
vertex skeleton diagram (which vanishes in the Abelian limit )|. Each of these three 
terms contributes starting at order a?. In this case, the skeleton decomposition of 
(|45| ) appears to be too naive: each of these skeleton terms has its own momentum 
fiow. In particular, if the BLM series is to have a one-to-one correspondence with the 
skeleton expansion one should write. 



ai?(<5^) = + sla{klp) + s1^a{kl.^p) + s\^a{kl^g) + 



(4^ 



To arrive at such a BLM series one should further decompose the coefficients in 
as follows. 



a + 



+ 



S2 + U 



Si + r^i^Po 



V2I + s 



, „3g (1) 



(49) 



where 



■'I 



(t'iCf 



a 



np 



(50) 



and 



2 



+ 



)Ca. (51) 



* Note that the three gluon vertex, which is a fundamental vertex in the theory, cannot be 
considered as just renormalizing the gluon propagator and the quark vertex. Part of it must define 
a new skeleton. This is in contrast to other diagrams appearing at this order which just renormalize 
the propagators or the quark vertex, and are therefore not candidates for new skeleton structures. 
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Here the combination [Cp — ^Caj corresponding to the non-planar skeleton (np) is 
suppressed in the large Nc liniitQ. Since the two Abelian parts of si, namely and s"^, 
are separately calculable, the coefficients r^^, r^Jip and r^Jg are uniquely determined 
from r2. Thus, in this example the color group structure plus the Abelian skeleton 
decomposition allow one to determine the non-Abelian skeleton decomposition. In 
the general case, where more skeleton structures are possible, this information will 
not suffice, and the decomposition of the coefficients will require a more complete 
understanding of the non-Abelian skeleton expansion. 

To summarize, we have seen that by tracing the flavor dependence of the pertur- 
bative coefficients in the skeleton scheme, one can identify the contribution of the 
different skeleton terms. This procedure allows us to "reconstruct" the skeleton ex- 
pansion algebraically from the calculated coefficients as summarized by eq. (^). This 
decomposition implies a unique BLM scale-setting which has a one-to-one correspon- 
dence with the skeleton expansion. We also learned that there are several limitations 
to the algebraic procedure which can probably be resolved only by explicit diagram- 
matic identification of the skeleton structures and the skeleton effective charge. These 
limitations include the need to 

a) treat separately contributions from skeleton structures which involve fermion 
loops (in the Abelian case these are just the light-by-light type diagrams) 

b) identify separately the different moments ^j-"'''^'* of a given momentum distribution 

function which appear as a sum (with any j and k such that j + k = n) in the 
perturbative coefficients of (3^ a^^*"^" 

c) identify separately the contributions of different skeleton terms which happen to 

appear at the same order in a. 



4 Using the ECH method in the framework of the 
skeleton expansion 

As we saw in the previous section, the essential ingredient of the BLM approach, 
which crucially relies on the skeleton expansion, is to disentangle running-coupling 
effects and treat them separately from the remaining expansion. Technically, this is 
realized by performing a skeleton decomposition. Running coupling effects can then 
be resummed in various ways, aiming at the approximation of the skeleton integrals 
Ri. First, there is the possibility to perform a full all-order resummation by eval- 
uating the integrals using some regularization in the infrared region (see e.g. |pT| ). 
This, however, requires the computation of the corresponding momentum distribution 
function, which can currently be done only at the level of the leading skeleton term. 

^In SU(A''c) the combination (Cf — jC'a) is sub-leading in Nc compared to Ca — Nc and Cp — 
{Nc'' ~ l)/{2Nc). 
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At the more modest level, the skeleton decomposition itself (|^) provides some infor- 
mation about the first few moments of the momentum distribution functions which 
is then used (section 3) to perform BLM scale-setting. Alternatively, the same infor- 
mation can be used to approximate the skeleton terms Ri{Q^) in the ECH method, 
which is particularly fit to deal with running coupling effects (see also [0). We 



shall see that this method has close relations with the scale setting procedure, but it 
also has some advantages over the latter. 

In this section we demonstrate how the ECH method can be used to provide 
resummation of running coupling effects in the framework of the skeleton expan- 
sion. The basic idea is that each skeleton term i?j(Q^) in our ansatz (H) is a 
renormalization-group invariant effective-charge raised to some power, i.e. one writes 
i?i(Q2) = (aR^(Q2))*+i instead of (||). Thus an^Q^) = (i?i(Q2))i/(^+i) ^an be simply 
be evaluated in the ECH method [Q, avoiding any explicit scale-setting procedure. 
In this method a/j. is computed by inverting the integrated renormalization-group 
equation, 

lnQVA% = (52) 
Jo /JR^ia) 

Finally, the observable aniQ'^) will be written as (cf. eq. (|^) and 



aR{Q^) = aR,{Q^) + siaR,{Q^f + S2aR,{Q^f + ■ ■ ■ . (53) 

Consider first the effective charge defined by the leading skeleton term = Rq 
as expanded in (^. From the next-to-leading order coefficient in this equation it 
follows that the ratio between the two scale parameters characterizing aR^^ and a is: 

AKo/A' = e-^^''. (54) 

This ratio is fully determined by the center of the momentum distribution function 
(which is also the leading order BLM scale-shift to.O; cf. eq. ([T^)) and is not mod- 
ified at higher orders. The latter affect just the corresponding ECH /3 function, 
Piioi^Ro) = dttR^/ InQ'^. Using the next-to-next-to- leading order expansion of qr^^ in 
terms of d and applying the general relation between effective charges 0, we have 

P^=P2 + Po{r2-rl)-P,n, (55) 

where /32 and are the three-loop /3 function coefficients of the skeleton coupling 
and of qr^, respectively. Using now ri and r2 of eq. (W) we obtain 



Pi = P2 + M^Pl (56) 



This means that for any momentum distribution 0o, (^2° is simply a sum of a universal 
piece /32, which characterizes the skeleton coupling, and an observable-dependent 
piece, namely the width of (po (see section 3) multiplied by /3q. 
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Recall that the three-loop j3 function coefficient in the skeleton scheme (32 is a 
polynomial of order 2 in /3o, namely = /^2,o + P2,iPo + (^2,2^0 (^^^ footnote 
following eq. @). Therefore is given by 



/3f 



P2,0 + P2,lp0 + P2,2p0 + 



^(2) 
' 2 



(57) 



In the large (3o limit is dominated by the last term, namely by the width of 
momentum distribution (pQ. In this case it is therefore the width which controls 
the convergence of the ECH j3 function, i.e. the accuracy of the calculated effective 
charge. Note that the same parameter controls the accuracy of the leading order 
BLM approximation [ll4| , |l^. Away from the large (3o limit, a small width implies 
proximity of and (32 (see Appendix A). Thus only if the universal (32 is not large, 
a small width implies smallness of i.e. good convergence of the effective charge 
approach applied to Ro- Similarly at the four-loop level, one gets 



^3 + 2Ms(3^ + 5M2(3i(3l 



{5t 



As usual 1^ the effective-charge a^jp is characterized by the scale ratio A|,^/A^ and 
(3 function coefficients. The same holds for higher skeleton terms. For example, it 
follows from (|36|) that aj^^ = (i?i)^^^ is characterized by A^j^A^ = e~5''2 and = 



P2 + 



,(2) 



(3q. Note that these are the same combinations appearing in 



the BLM scale-shift for Ri, eq. (|4|). 

It is also worth noting that the suggested effective charge approach yields a result 
identical to the BLM scale setting method applied in the skeleton scheme, in the 
approximation where the (3 functions of the effective charges associated to the various 
skeletons (3^^ are all replaced by the skeleton coupling (3 function, (3. This is equivalent 
to assuming that, except the average, all the central moments of the characteristic 
functions vanish identically. Then both approaches effectively yield a multi-scale 
series where the scales correspond to the average momentum flowing in each skeleton 
diagram, kifl. 

To conclude, we have shown that the explicit scale-setting procedure can be re- 
placed by the ECH method. One advantage is that the latter does not suffer from 
the scheme and scale ambiguities still present (see the footnote before eq. (|20|)) in 
the series for the BLM scales. We stress that our example here heavily relies on the 
specific ansatz assumed for the skeleton expansion. However, contrary to the BLM 
scale setting method, the suggested effective charge approach would apply equally 
well to more general cases where e.g. the two couplings in eq. (0) are different. 



5 BLM and conformal relations 

Let us now consider the general BLM scale-setting method, where the scheme is 
not necessarily the one of the skeleton effective charge, and no correspondence with 
the skeleton expansion is sought for. Then any scale-setting procedure which yields 
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an expansion of the form (|T^) with Nf independent q coefficients and scale-shifts 
which are power series in the couphng is legitimate. We saw that under these 
requirements there is no unique procedure for setting the BLM scale beyond the 
leading order (fco,o)- Nevertheless, as we now show, the coefficients Cj are uniquely 
defined. In fact, the q have a precise physical interpretation as the "conformal 
coefficients" relating and a in a conformal theory defined by 

/3(a) = -Poa^ - pia^ + ■ • ■ = 0. (59) 

To go from real-world QCD to a situation where such a conformal theory exists 
one has to tune Nf-. when Nf is set large enough (but still below ^^c, the point 
where asymptotic freedom is lost) (3i is negative while jSo is positive and small. Then 
the perturbative (3 function has a zero at app ^ — /5o//3i; i-e. there is a non-trivial 



infrared fixed-point The perturbative analysis is justified if /5o, and hence 

Opp, is small enough. 

Physically, the existence of an infrared fixed-point in QCD means that correla- 
tion functions are scale invariant at large distances. This contradicts confinement 
which requires a characteristic distance scale. In particular, when /3o — > the in- 
frared coupling is vanishingly small. Then it is quite clear that a non-perturbative 
phenomenon such as confinement will not persist. The phase of the theory where the 
infrared physics is controlled by a fixed-point is called the conformal window. In this 
work we are not concerned with the physics in the conformal window^. We shall just 
use formal expansions which have a particular meaning in this phase. 

The BLM coefficients q are by definition A^j-independent. Therefore the expan- 
sion of aji according to eq. (p!8| ) is valid, with the same Cj's both in the real world 
QCD and in the conformal window. In the conformal window a generic coupling 
a(A;^) fiows in the infrared to a well-defined limit a(A;^ = 0) = Opp. In particular, 
eq. ([ISD becomes 

o-r" = %p + Ciflpp^ + C2app^ + CsOpp'^ + ■■ ■ (60) 

where we used the fact that the /cj's are proportional to Q, which follows from their 
definition = Q'^exp{—ti), together with the observation that the scale-shifts in 
( [23| ) at any finite order are just constants when a^Q"^) — > a^p. Eq. (^) is simply the 
perturbative relation between the fixed-point values of the two couplings (or effective 
charges) and a. 

Note that in this discussion we ignored the complication discussed at the end 
of section 3, concerning the possibility of applying BLM scale-setting in the case 
of several skeletons contributing at the same order (cf. eq. (0)). In this case the 
argument above holds as well, while the conformal coefficients will be the sum of all 
BLM coefficients appearing at the corresponding order. For the example considered 
in section 3, (|49|) , we would then have ci = si = + s"^ + sl^ . 

this phase is investigated from the point of view of perturbation theory in both QCD 
and supersymmetric QCD. 
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According to the general argument above, the BLM coefficients (|T8[) should coin- 
cide with the conformal coefficients in dBD]). In the next section we calculate conformal 
coefficients directly and check this statement explicitly in the ffist few orders. 



6 Calculating conformal coefficients 

Let us now investigate the relation between the conformal coefficients q appearing 
in (pOD and the perturbative coefficients r^. 

For this purpose, it is useful to recall the Banks-Zaks expansion: solving the 
equation /3(a) = in (|59| ) for such Nf where Pq is small and positive and Pi is 
negative, we obtain: a^p ^ —Po/Pi > 0. If we now tune Nf towards the limit 
Y^c from below, (3q and therefore Cpp become vanishingly small, which justifies the 
perturbative analysis In particular, it justifies neglecting higher orders in the 

f5 function as a ffist approximation. In order to take into account the higher orders 
in the (3 function, one can construct a power expansion solution of the equation 
P{a) = 0, with the expansion parameter as the leading order solution, 

«o = --r, = —5-- (61) 

Pl|/3o=0 ~/^l,0 

In the last equality we defined /?i = Pi^ + Pi^iPo where Pij are A^/-independent. 
Similarly, we define^ for later use 

P2 = P2,0 + P2,lPo + P2,2Po + P2,3Pl (62) 

We shall assume that the coupling a has the following Banks-Zaks expansion 

Opp = ao + Vial + V2al + fsag H (63) 

where Vi depend on the coefficients of P{a), see e.g. [^]. For instance, the ffist 
Banks-Zaks coefficient is 

vi = - (64) 

Pi,o 

Suppose that the perturbative expansion of aniQ'^) in terms of a(Q^) is given by 
aniQ^) = a{Q^) + naiQ^f + r^aiQ^f + ■■■ (65) 

Based on the fact that are polynomials of order i in Nf, and that ao is linear 
in Nf, one can uniquely write a decomposition of rj into polynomials in oq with 
iV/-independent coefficients 

ri = r-1,0 + ri,iao (66) 
r2 = r2fi + r2,iao + r2,2ao^ 



§We recall that in the skeleton scheme P2.3 — 0. 
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and so on. For convenience we expand here in oq rather than in (3o. The relations 
with the "skeleton decomposition" of ri and r2 in eqs. (|T^) and (PT| ) (or in are 
the following 

ri,i = -/^i^or^) r2,i = -PiA'^r^^^ - Pi,oPiA^ (67) 



r2,2 = Plor; 



For ra we have, based on 



2 



r-s.o = rf + rJ'Vf ^A.o + ^/^a.o- (68) 



Using eq. ( |65D at = with (^) and the Banks-Zaks expansion for Cpp (§^), it 
is straightforward to obtain the Banks-Zaks expansion for a^^ 

an" = «o + ujio-l + W2al + Wsa^ H (69) 

with 

wi = vi + ri^o (70) 
W2 = V2 + 2rifiVi + ri,i + r2,o 

W3 = V3 + 2ri^oV2 + Ti^oVi"^ + 2ri,it;i + 3r2,o^^i + ^2,1 + rg^o 

Having the two Banks-Zaks expansions, one can also construct the series which 
relates two effective charges aff and Opp at the fixed-point. Inverting the series in 
([63| ) one obtains Qq as a power series in Opp, 

ao = Opp + MiOpp + M2app + -usflpp + ■ ■ ■ (71) 

with Ml = — t>i and U2 = vf — f 2 etc. Substituting eq. (^) in ( |5^ one obtains the 
"conformal expansion" of a]^ in terms of app according to eq. (^) with 

ci = ri,o (72) 

C2 = + r2,o 

C3 = -^^1,1^ 1 + r2,i + r3,o 

C4 = 2ri,ii;i^ - ri,it;2 - r2,iVi + r2,2 + rs^i + r^^ 



Thus the coefficients Vi of the Banks-Zaks expansion ( |63D and the coefficients rj of 
(|65| ) are sufficient to determine the conformal coefficients q to any given order. 

Clearly, the Banks-Zaks expansions (|63|) and (|69D and the conformal expansion of 
one fixed-point in terms of another (^) are closely related. Strictly speaking, both 
type of expansions are meaningful only in the conformal window. However, we saw 
that the coefficients of (|60D coincide with the ones of the BLM series ( plSf ) which is 
useful in real world QCD. We recall that the general argument in the previous section 
does not depend on the specific BLM scale-setting prescription used, provided that 
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the scales ki are proportional to Q and the q's are Nf independent. Comparing 
explicitly Ci, C2 and C3 in eq. (|72|) with the BLM coefficients obtained in the previous 
section, namely Cj = r|°\ we indeed find that they are equal (compare using eq. (p7|), 
)). In particular, the "skeleton decomposition" of eq. (^51) , which can be 
formally performed in any scheme, provides an alternative way to compute conformal 
coefficients. 



m) and (64 



7 Examples 

The skeleton expansion assumption implies that the skeleton (conformal) coefficients 
Si are free of running coupling effects. In particular, contrary to the standard per- 
turbative coefficients in a standard scheme such as MS, the large order behavior of 
conformal coefficients is not dictated by renormalon factorial increase, and should 
therefore be softer. 

In other words, the effective convergence of the fixed-point relation (^) where 
a is taken as the skeleton coupling effective charge a is expected to be better than 
standard perturbative expansions. As we shall see in section 7.3, this expectation is 
not restricted to the skeleton scheme but applies also to general conformal relations, 
e.g. between two physical effective charges. In addition, if we assume that the 
skeleton coupling (3 function itself is renormalon-free, it follows that also the Banks- 
Zaks expansion of a generic physical quantity an is renormalon-free. This is because 
the latter assumption implies that the Banks-Zaks expansion of d (^) is free of 
renormalons, and then, by substituting it in the renormalon-free conformal relation 
between the observable and a, one recovers the Banks-Zaks expansion of a/?, which 
must therefore be renormalon-free as well. 

Thus, the general expectation is that all conformal and Banks-Zaks relations 
are free of renormalons and have better convergence properties. Our purpose here 
is to examine through available examples in QCD whether this expectation is re- 
alized. Indeed, as we recall below, it has been noted by several authors (e.g. in 



12| , |33| , p!q , p6[ |37| ) that conformal coefficients and Banks-Zaks coefficients are typ- 
ically small. We would like to interpret these observations based on the assumed 
skeleton expansion and relate them to the absence of renormalons. As concrete ex- 
amples we shall concentrate on the following observables: 

a) The Adler D-function, 



^2rfn(Q2)_ 3 



1 + - Cpao 
4 



(73) 



where a^j is normalized as an effective charge, and ^{Q^) is the electromagnetic 
vacuum polarization, 

Ai^^i J rf^xe^'?-^ (0|T{j'^(x),j^(0)} |0) = (gV - ^{Q^)- (74) 
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b) The polarized Bjorken sum-rule for electron nucleon deep-inelastic scattering, 



dx 



9A 

6 



(75) 



c) The non-polarized Bjorken sum-rule for neutrino nucleon deep-inelastic scatter- 
ing, 



dx 



Frix,Q')-Frix,Q' 



1 

= 1 -—ap,. 



d) The static potential, 



ay 



(76) 



(77) 



In all four cases perturbative calculations have been performed (refs. [^] through 
PB| , respectively) up to the next-to- next-to-leading order r2 in eq. (|T]). 

For later comparison with conformal relations, we quote some numerical values 
of the coefficients in the standard perturbative expansion in Gms = Oms (Q^) for the 
vacuum polarization D-function (|73D 



(7^ 



O'D — ^Ms + di a|jg + 
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2.0 


18.2 




Nf = 
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1.6 


6.4 




Nf = 


= 3 


0.14 


-27.1 




Nf = 


= 16 


1.06 


14.0 




Nf- 


= 0..16 


the polarized Bjorken sum-rule (^ 


I) 










k2 


n— 4- ■ 






4.6 


41.4 




Nf 


= 


3.5 


20.2 




Nf 


= 3 


-0.75 


-34.8 




Nf 


= 16 


2.1 


21.0 




Nf 


= 0..16 



(79) 



where in the first three lines in (|7q) and (^) the coefficients are evaluated at given 
Nf values, while the last line corresponds to an average of |rj| in the range Nf = 
through 16. 

We see that the coefficients in a running coupling expansion in the MS scheme 
increase fast already at the available next-to-next-to-leading order. This increase has 
been discussed in connection with renormalons, for example in @]. A priori, it is hard 
to expect that the large-order behavior of the series will show up already in the first 



few leading orders. We mention, however, that in ref. the Bjorken sum rule series 
(for Nf = 3) was analyzed in the Borel plane based on the three known coefficients, 
indicating that the first infrared renormalon at p = 1 does show up. 
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7.1 The Banks- Zaks expansion 

Let us now compare the magnitude of the coefficients in the standard expansion, e.g. 
in eqs. (^) and ([79|), to that of conformal coefficients. For the latter, one can choose 
to examine conformal relations between effective charges (see section 7.3) or - the 
Banks-Zaks expansion. 

The Banks-Zaks expansion for the fixed-point value of the vacuum polarization 
D-function ( |75D is 

af^ = ao + 1.22 ao^ + 0.23 ao^ + ■•• (80) 

whereas for the Bjorken sum-rule it is 

a^p= ao + 0.22 ao^ - 1.21 ao^ + (81) 

Comparing (^) and (|80| ) with the corresponding running coupling expansions in 
MS, namely (|78D and (|79D, the difference in magnitude of the coefficients is quite 
remarkable [0, ^ . Taking into account the fact that the coefficient of a}^^ contains. 



among other terms, a Ca' = 3* term, this fast apparent convergence seems rather 
surprising. From this point of view, the absence of renormalons may not be considered 
a sufficient explanation. 

For the non-polarized Bjorken sum-rule defined by (|82D, the Banks-Zaks coeffi- 
cients are even smaller 

aF;p = ao - 0.45 ao^ + 0.16 oq^ + ••• (82) 

and exhibit an impressive cancelation of numerical terms appearing in the running 
coupling coefficients |^ . The static potential shows a different behavior. In this case 
the Banks-Zaks expansion [R^, R3] 



a^p = ao - 0.86 ao^ + 10.99 ao^ + ■■• (83) 

has a significantly larger next-to- next-to-leading order coefficient. Taking into ac- 
count the numerically large color group factor Ca^ = 9, the magnitude of this next- 
to- next-to-leading order coefficient is quite reasonable. 

Another physical quantity for which the Banks-Zaks coefficients are relatively 
large is the critical exponent 7 [0, |3^, |3^ 



Po da 
where 



(84) 

a=app 



7= ao + 4.75 ao^ - 8.89 oq^ + ••• (85) 

Since this quantity does not depend on Q^, there is no direct comparison between a 
running coupling expansion and the Banks-Zaks expansion. 

To conclude, we have seen that the Banks-Zaks coefficients for physical quantities 
typically have smaller coefficients compared to the standard running coupling expan- 
sion. In some cases, their convergence is surprisingly good, even taking into account 
the absence of running- coupling effects. 
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7.2 Conformal relations in the skeleton scheme 



Examining the Banks-Zaks expansion we found that the coefficients are significantly 
smaller than standard running-coupling coefficients. The same conclusion would fol- 
low from examining direct conformal relations between observables. The coefficients 
of such relations (see section 7.3) are not only small, but also exhibit a remarkable 
simplicity [Q. Both the smallness and the simplicity of these coefficients seem a 



natural consequence of the conformal limit. The smallness, in particular, is naturally 
attributed to the absence of running-coupling effects. 

The first step in trying to substantiate this statement in the framework of the 
postulated skeleton expansion is to consider the conformal relations in the skeleton 
scheme. It is natural to expect that conformal relations between observables and the 
skeleton coupling will be small, thus explaining the above observations. 

To this end, let us consider now the conformal relation in the skeleton scheme (l60|) 
as defined by the pinch technique. Since the skeleton coupling a has been identified 
only at the one-loop level (p!2D , our information on the coefficients Si is quite limited: 
by a direct calculation (using the next-to-leading order coefficient ri and either ( [T^) 
or (0)) we can only determine Si. For example, for the observables defined above it 
is 

r -(l/4)C^ - {1/S)Cf = -11/12 D 
, _ JO) _ -(1/4)C^ - (7/8)C^ = -23/12 g, 

~ 1 ~ I -(1/4)C^ - (11/8)C^ = -31/12 Fi ^^"^ 
[ -Ca = -3 V 

Note the absence of a Ci;' term in the case of the static potential. This can be 
understood based on the Abelian limit, where it is known that this effective charge 
coincides with the skeleton coupling (there, the Gell-Mann Low effective charge) up 
to light-by-light type corrections. Therefore the momentum distribution function of 
the leading skeleton term 0o is just a (5-function, 0o(^^) = and in the Abelian 

limit there are strictly no (A^j-independent) sub-leading skeleton terms. 

The higher-order coefficients Sj, for i >2, depend on yet unknown characteristics 
of the skeleton coupling scheme. In particular, as we discuss in Appendix A, S2 



depends on the skeleton (3 function coefficient (32- However, as can be seen in eq. ( |105|) 
there, the dependence on this coefficient cancels in the difference of S2 between any 
two observables, which is therefore calculable. 

Without a diagrammatic identification of the skeleton structure, one cannot isolate 
skeletons with fermion loops attached to three gluons, which may appear at the order 
considered. Therefore we shall just treat the entire Nf dependence (excluding Abelian 
light-by-light diagrams) as if it appears due to the running coupling, according to 
eq. (145|) where S2 is Nf independent. For the observables considered above we then 



find: 



3 ^ ^ 3 
4 
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■ vr 



Ca' 



7.045 

25 



(87) 



16 



23 
32 



19.66 



This gives some estimate of the size of S2 for these observables. The S2 coefficients 
turn out to be larger than the Banks-Zaks coefficients quoted above (as well as the 
conformal coefficients in the relation between observables). They can even be compa- 
rable in size to the next-to-next-to- leading order coefficients in MS. Thus the assumed 
form of the skeleton expansion does not provide a satisfactory explanation for the 
observed smallness of conformal coefficients. We stress again that we mistreated here 
the Nf dependence which is associated with the skeleton structure, namely fermion 
loops attached to three gluons. Eventually, this will have some impact on the mag- 
nitude of the (A'^^ dependent) skeleton coefficients S2, which we cannot evaluate at 
present. 



7.3 Direct relations between observables 

As we saw above the knowledge about the skeleton coefficients is very limited beyond 
next-to-leading order. However, there is a way to consider systematically conformal 
relations avoiding the use of the skeleton scheme. Having renormalon-free conformal 
expansions (|60|) for two QCD observables in terms of the skeleton effective charge a, 
one can eliminate the latter to obtain a direct conformal relation between the two 
observables. The existence of a skeleton expansion (Q) for the two observables implies 
that this conformal relation is free of renormalons. 

Conformal coefficients of this type can be computed either from the Banks-Zaks 
expansion (|7^ ) or in the framework of BLM, as the coefficients in a commensurate 
scale relation |]T2[. The latter can be obtained by applying BLM directly to the per- 
turbative relation between two observable effective charges (and so it does not require 
identification of the skeleton coupling). However, it should be noted that whereas 
the one-to-one correspondence between the BLM series and the skeleton expansion 
specifies a unique scale-setting procedure when the skeleton scheme in used, the 
scale-setting procedure in direct relations between observables remains ambiguous. 
As explained in sections 3 and 4, the conformal coefficients themselves are uniquely 
determined, independently of the particular way the scales are set. 

In addition to being numerically small, conformal coefficients in the direct rela- 
tions between observables turn of to be simpler |1^, in terms of color group factors 



and numerical (n terms. This simplicity is naturally attributed to the conformal 
limit. 

There is one example where a direct all-order conformal relation is known - this 
is the Crewther relation relating the vacuum polarization D-function effective charge 
ao, defined by (|75D, with the polarized Bjorken sum-rule effective charge a^^, defined 
by (^). The Crewther relation is 



3 

a-gi - flD + -^Cpag^aD = -(3{a)T{a) 
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(89) 



where T{a) is a power series in the couphng 

T(a) = Ti + - 

and Ti are polynomials in Nf. 

If has a perturbative fixed-point a^jj, then it is convenient to write the 
r.h.s. of ( pSD in terms of a^j. /?(a^) = and so the r.h.s. vanishes at = 
corresponding to the infrared limit. Therefore a^j also freezes perturbatively, leading 
to the original conformal Crewther relation 



''D 



91 



1 + f Cf^d' 



(90) 



Taking Nc = 3 we have Cp = | and then the conformal coefficients are just one to 
any order in perturbation theory, 



91 



+ {<iT+{<j 



(91) 



Being a geometrical series this conformal relation provides a nice example of a pertur- 
bative relation free of renormalon divergence. In addition, it exemplifies the simplicity 
of the conformal limit: here the conformal coefficients do not contain any non-Abelian 
Ca terms. 

As noted in (see also [Q) it is possible to write for two generic observables 
A and B, at two arbitrary scales Qa and Qb, the following decomposition of the 
perturbative series relating the two. 



aA = CABiaB) + Pia'B)TAB{aB) 
Here Cab is the "conformal part" of the series, i.e. 

CAsiaB) = aB + cia| + C2a| + ■ 



(92) 



(93) 



where q are the conformal coefficients appearing in the expansion of in terms 
of a^, and Tab{o,b) is a perturbative series of the form (|89|). In other words the 
non- conformal part of the relation between the two observables is factorized 
as (3{aB)TAB{o.B)- Taking the limit /5 — > then gives the conformal relation. In 
particular, one can write such a factorized relation between an observable effective 
charge and the skeleton coupling. Then the conformal coefficients q in (p3|) are 
the skeleton coefficients Sj. Explicitly, this can be shown based on the skeleton 
decomposition of the series 



an 



a + Sio^ + 320? + s^a'^ + ■ ■ 



+ 



(94) 



X 



Finally, we also quote the conformal relations between the vacuum polarization 
D-function and the non-polarized Bjorken sum-rule (|7BD, 

,2 . / „„\3 



+ 1.67 (a^^) +1.57(a^^) +■ 



(95) 
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as well as the static potential ([77|) , 

= + 2.08 (a^^)' - 7.16 (a'y^f + ■■■. (96) 

Taking into account the C^* = 3* contribution to Cj, these expansions all seem well- 
behaved. 



7.4 Expansions in MS 

Finally, it is interesting to return to the expansion in MS and examine the corre- 
sponding conformal relations. Such relations turn out to have large coefficients. For 
example, 

a- = - 0.083 (a^ f - 23.22 (a^ )=^ + ■ ■ • (97) 

and 

< = «^ - 0.917 (a^ f - 22.39 (a^ )=^ + ■ ■ • (98) 

have large next-to-next-to-leading order coefficients, in a striking contrast with the 
conformal relation (0) between a^J and a^^. Note that these large conformal coeffi- 
cients do not provide an explanation of the large coefficients in ( fTSD and (|79|). The 
former are by assumption independent of Nf, as opposed to the latter. For small 
Po (e.g. Nf = 16) the negative sign (and eventually also the magnitude) of the full 
coefficient can presumably be attributed to the conformal part. However, for larger 
values of /3o, relevant to real world QCD, the non-conformal part clearly dominates 
making the full next-to-next-to-leading order coefficients positive. 

These large conformal coefficients in (^^ and (^8]) are due to an intrinsic property 
of the MS coupling, since they appear already at the level of the Banks-Zaks expansion 



37, 38 



= ao + 1.1366 ao^ + 23.2656 oq^ + (99) 

Note that has, by far, a larger next-to-next-to-leading order Banks-Zaks coeffi- 
cient compared to any known physical effective charge. 

We stress that the large next-to-next-to- leading order coefficients in (p?]), ( p^ ) 
and (p9|) are not associated with renormalons. The MS f3 function, being defined 
through an ultraviolet regularization procedure, should not be sensitive to the in- 
frared. Therefore infrared renormalons are not expected. It is more difficult to draw 
any &m conclusion concerning the absence of ultraviolet renormalons. Since there 
seems to be no reason to assume a skeleton structure or any other representation in 
the form of an integral over a running coupling, we suspect that ultraviolet renor- 
malons do not exist there as well. 

To conclude, the case of conformal relations in MS teaches us not to associate 
automatically any large coefficient in QCD with running-coupling effects. Indeed, 
in field theory there are other sources of large coefficients, such as multiplicity of 
diagrams. 
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8 Conclusions 



The fast growth of perturbative coefficients and the related renormahzation scale 
and scheme ambiguities of perturbative expansions have greatly limited the predic- 
tive power of QCD. In many cases, this divergent behavior is predominantly due 
to running-coupling effects. The existence of an Abelian-like skeleton expansion in 
QCD would make it possible to disentangle in a unique way such effects, separating 
them from the conformal part of the perturbative expansion of a generic physical 
quantity. The effect of the running coupling could then be treated systematically to 
all orders in perturbation theory in a renormalization-scheme invariant manner by 
renormalon-type integrals. The normalization of these skeleton integrals is controlled 
by conformal coefficients that are hopefully better behaved, making the truncated 
skeleton expansion a better approximation to the physical observable compared to 
the standard perturbative expansion of the same order. 

Resummation of running coupling effects has in many cases a significant role in 
phenomenology . Direct resummation is currently restricted to the level of a single 
dressed gluon, where the Abelian large Nj limit can be used. The formulation of 
perturbation theory in the form of a skeleton expansion has implications which go 
beyond the perturbative level. In particular, it provides a natural framework to deal 
together with the resummation and the related power-corrections. The renormalon 
integral contains essential information on the type of power-corrections one should 



expect for a given observable. Moreover, it can be used to combine ||3y, |31| such power 
corrections with the perturbative expansion avoiding double-counting or dependence 
on the particular prescription used to regularize infrared renormalons. These aspects 
were discussed in detail in for the example of the average thrust. 

In this paper we have concentrated on the conformal part of the perturbative 
expansion, based on a postulated ansatz for the skeleton expansion. We have shown 
that the (A^j-independent) coefficients of this expansion and of the related BLM se- 
ries have a precise interpretation when a perturbative infrared fixed-point is present: 
they are the conformal coefficients in the series relating the fixed-point value of the 
observable under consideration with that of the skeleton effective charge. The per- 
turbative infrared fixed-point appearing in multi-fiavor QCD allows one to calculate 
these conformal coefficients through the Banks-Zaks expansion. We stress that the 
identification of the skeleton coefficients with the ones of the conformal relations de- 
fined in the small /Sq limit strongly relies on the particular ansatz we have taken, 
namely that the entire Nf dependence originates in the running- coupling itself, leav- 
ing the conformal coefficients Nf independent. On the other hand, the identification 
of the BLM coefficients with those of the conformal relations of the small (3o limit does 
not rely on any additional assumption, and it holds independently of the particular 
way BLM scale-setting is performed. 

Existence of an underlying skeleton structure implies that BLM (conformal) co- 
efficients do not diverge factorially due to renormalons. Of course, there can be 
other effects which could make these coefficients diverge such as combinatorial fac- 
tors related to the multiplicity of diagrams. Since in QCD this type of divergence 
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is much softer than that of renormalons, we expect the BLM and possibly also the 
Banks-Zaks expansions to be "better behaved". This expectation is supported to 
some extent by previous observations concerning the smallness of the first few known 



BLM coefficients I1T2| and the Banks-Zaks coefficients [03, WA |38|. On the other hand 



the absence of renormalons does not always seem to be a sufficient explanation of 
the observed difference between conformal and non-conformal coefficients. At the 
same time, large coefficients which are not associated with running- coupling effects 
do appear in QCD, e.g. in conformal relations with the MS coupling. 

The uniqueness of the skeleton coupling in QED, which is identified as the Gell- 
Mann Low effective charge, is an essential ingredient of the dressed skeleton expan- 
sion. It is still an open question whether an Abelian-like skeleton expansion exists 
in QCD and what the constraints are which would determine the skeleton coupling 
uniquely. The pinch technique may provide the answer [El, 121, |2H| once it is sys- 



tematically carried out to higher orders. We recall that the skeleton coupling is not 
constrained from the considerations raised in this paper: the only requirement fol- 
lowing from the large Nf limit is that Pi in this scheme does not contain an Nj'^^ 
term. Since the decomposition of the coefficients (|45|) can be performed in any scheme 

(i) 

yielding the moments to arbitrary high order, the corresponding functions 0i can 
be formally constructed, up to the limitations discussed in section 3.4. It thus seems 
that one can formally associate a "skeleton expansion" to any given coupling. The 
absence of renormalons in the conformal coefficients in a specific scheme implies that 
there are other schemes which share the same property: it is straightforward to see 
from the definition of the skeleton terms Ri that an A^j-independent re-scaling of the 
argument of the coupling leaves the conformal coefficients unchanged. More gener- 
ally, any "renormalon-free" transformation of the skeleton coupling would leave the 
"skeleton coefficients" free of renormalons. It is certainly interesting to find further 
constraints on the identity of the skeleton effective charge in QCD. 

The BLM method provides a pragmatic way to deal with running-coupling effects 
beyond the single dressed gluon level. By decomposing the perturbative coefficients in 
the specific way implied by the skeleton expansion the contributions from the different 
skeleton integrals as well as the conformal coefficients can be identified. The BLM 
scales are then set such that there is a one-to-one correspondence between the terms 
in the BLM series and the skeletons, provided that BLM scale setting is performed 
in the skeleton scheme. As an alternative to the BLM scale setting procedure we 
saw that the skeleton integrals can also be approximated by applying the method of 
effective charges to the separate skeleton terms. 

In practice, BLM scale-setting can also be applied in physical schemes yielding 
a commensurate scale relation. This way conformal relations, which have a natural, 
maximally convergent, form (like the conformal Crewther relation) can be used as a 
template for real-world QCD predictions, even if the underlying skeleton structure is 
not completely understood. Still, since the conformal relation between the fixed-point 
value of a generic observable and that of the skeleton effective charge is renormalon- 
free, it follows, upon eliminating the skeleton effective charge, that the coefficients in 
commensurate scale relations between observables are also renormalon-free. When 
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such a conformal template is used for real-world QCD calculations, the effect of the 
non-zero /3— function is to modify the values of the scales k"^ of the effective charge at 
each order of the expansion. We stress, however, that having no more correspondence 
with the skeleton expansion, the motivation for a particular scale-setting procedure 
is lost. 

The BLM procedure cannot replace an eventual diagrammatic formulation of the 
skeleton expansion. We saw that the scale-setting prescription depends on the ansatz 
for the skeleton expansion, and any unknown concerning the form of the latter would 
have some impact on the former. We have considered several ways in which the 
simple ansatz we introduced (0) may be generalized. This includes in particular the 
possibility that several skeleton diagrams will appear at the same order and that 
certain skeleton diagrams will contain some fermion loops as part of their structure, 
making the corresponding conformal coefficients Nf dependent. In addition, we have 
seen that even in the case of a simple form of the skeleton expansion, the skeleton 
decomposition of the coefficients, cannot be performed up to arbitrarily high order 
just based on the Nf dependence, but rather requires some additional knowledge 
based on an explicit diagrammatic formulation. One should also be aware of the 
possibility that an Abelian-like skeleton expansion with a single effective charge might 
fail to exist in QCD. The non-Abelian skeleton expansion may then be based on 
several dressed Green functions, namely several different effective charges. Even in 
this more complicated case the most important properties of the skeleton expansion 
assumed here may hold. This includes the possibility to associate running-coupling 
effects to the various skeleton terms in a renormalization-group invariant way, and the 
interpretation of the skeleton coefficients as conformal coefficients when a perturbative 
infrared fixed-point is present. 
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Appendix A — The skeleton expansion and the ef- 
fective charge approach 

A priori, the skeleton expansion approach, which relies on the assumption of a uni- 
versal skeleton coupling, seems antagonist to the original effective charge approach 
[|5| which treats all effective charges independently and in a symmetric manner. In 
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section 4 we saw how the ECH method can be used in the framework of the skeleton 
expansion to approximate separately each skeleton term Here we revisit the orig- 
inal ECH approach which attempts to evaluate the entire observable directly, and 
examine it from the point of view of the assumed skeleton expansion. 

We begin by comparing the original ECH approach to the application of the 
ECH method for the leading skeleton term Rq. The first difference is, of course, 
in the ECH scale parameter. To facilitate the comparison, suppose that we start 
with a perturbative expansion (P of the observable in terms of a{Q^), with the 
corresponding coefficients rj. In the original ECH approach this implies a scale ratio 
of A|./A^ = e"''^/^". This can be compared with (|5^). The difference between the 
two is due to the rf^^ component in the next-to-leading coefficient ri, the component 
which is not associated with the leading-skeleton. In practice, in many cases in QCD 
the running-coupling component dominates the next-to- leading coefficient. In such 
cases the two scales are close. 

Next, also the f3 function of the ECH method, (3R{aR) = dan/lnQ'^ is different 
(beyond the universal two- loop order) from that of Rq. At the three-loop level the 
latter is given in (|56|) whereas the former is 



— (^2,0 + l32,iPo + (^2,2(^0 + 



J2) / (1) 
o '1 



Pi (100) 



where we exhibited the fact that the term leading in is the same in and 
■ As noted in section 4, in the large /3o limit is proportional to the width of 



the distribution (po, namely to 



- (r«^'' 



This remains correct also for the (3 



2 -yi J 

function of the full effective charge aji since adding sub-leading skeleton terms would 
not modify the leading O (/3q) term. For the four examples considered in section 7, 
this parameter is given in table I. 



/52,3 


/5?3 


/3f,3 









2.625 


2.389 


1.500 






Table 1: Comparison of effective charge (3 function coefficients in the large /3o ap- 
proximation given by the width of (f)Q, /32,3 = — (j'l^^^ ■ 

It is natural now to consider the possibility that Rq is a good approximation to 
the observable an. In the effective charge approach at the next-to-next-to-leading 
order, this can be realized if p^" is a good approximation to (32 ■ large (3q limit 

the two are equal. Beyond the large (3q limit one can ask whether 

/?2'!o+/52Vo+/52Vo ^ ^2=^2,0 + ^2,l/3o + ^2,2/3o', (101) 

namely whether (32 q + (32i(3q + (^2 2^1 ^ generic observable which admits a skeleton 
expansion is approximately universal and close to the three-loop skeleton coupling (3 
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function coefficient j32- If this holds for arbitrary f3o then 



2A 



(102) 



for i = 0, 1, 2. The violation of the equalities in (|101|) and (|102|) is, of course, due to 
sub- leading terms in the skeleton expansion Ri and i?2- This can be seen explicitly 
by substituting of eq. (|45|) in the general relation 



(5^='P2 + f^o{r2-rf)-P 



in 



to obtain the "skeleton decomposition" of (3. 



R 
2 ) 



= P2+ {s2 - sl) /3o + .1 (r« - 2r«) p', + [r? - (r«) 
Finally, decomposing P2 and Pi in terms of Pq, we obtain^] 



(103) 



/3o' - siPi. (104) 
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02,1 - Pl,lSi + (S2 - sfj Po 
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J2)_^i)y 



(105) 



Clearly, if for a given observable the skeleton coefficients determining the normaliza- 
tion of the sub-leading skeleton terms (sj) are small, then even away from the large 
Po hmit p^ will be close to 

In order to check ( |102| ) explicitly for a given observable, one needs to calculate the 
P function coefficients of both the observable effective charge and the skeleton 
effective charge P2,i. For the latter we currently know only /32,o (see below) and so 
the examination of ( [L02|) for P2,i and /32,2 cannot yet be accomplished. 

To obtain /?2,o we can use the general result |^ or, alternatively use eq. (|105|) , 
which is valid for a generic effective charge which admits a skeleton expansion. The 
latter yields, 

P2,0 = P2,0 + Pl,0Sl- (106) 

Using this relation for various effective charges, e.g. the vacuum polarization D- 
function ( [73D or the Bjorken sum-rule (ffSD, in the skeleton coupling scheme dlT 
defined through the pinch technique, we obtain 



Ca 
512 



(44^2 



88CylC p 



301CI 



and for Nr 



P. 



2,0 



26845 
1536 



-17.477. 



(107) 



108) 



^The scheme of the skeleton couphng can be parameterized at the three-loop order |5) by the 



next-to-leading order coefficient (si and rj^-*) and by (32 i.e. /92,i for i = 0, 1, 2. Eq. (|105|) then shows 
explicitly that the effective charge /3 function coefficient /3|^ determines uniquely the remaining coef- 



ficients of the "skeleton decomposition" ( [45D namely, S2, t^'^ and r2^''. This reflects the observation 
in section 3 that formally, the "skeleton decomposition" can be performed in any scheme. 
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Finally we check to what extent the suggested universality of the effective charge 
f3 function coefficients ( |102| ) holds for the four effective charges examined in section 6, 
namely the effective charges related to the vacuum polarization D-function (^) and 
the polarized ( [75| ) and non-polarized ( [76| ) Bjorken sum-rules, as well as the static 
potential. The known coefficients are listed in the following table. 



i 






/30 


m 


PI 





-17.477 


-23.607 


-30.294 


-34.753 


-37.54 


1 


? 


-16.032 


-11.282 


-6.903 


5.366 


2 


? 


8.210 


8.057 


8.783 


11.740 



Table 2: Comparison of effective charge /3 function coefficients. 

Although the coefficients i these observables have some common trend (e.g. 
for a given i the signs are the same, with the exception of for i = 1) it turns 
out that the fluctuations in their magnitude are rather large. In particular, in case 
of (32 which we know the value of the universal piece characterizing the skeleton 
coupling (32 fi, the latter and the contribution of the sub- leading skeleton Ri (through 
Si in eq. ( [105| )) are of the same order of magnitude. The fluctuations between different 
observables are moderate only for (32 2- 

In it has been observed that for the observables considered above (the static 
potential excluded) exhibit very close numerical proximity, especially for Nf = 



through 7. The extent to which universality of the sort examined here ( [102| ) holds is 
not enough to explain this flnding of ||3^ . 



The proximity of (32 2 ^^e various effective charges implies that applying multi- 
scale BLM scale-setting for one observable in terms of another, the second scale-shift 
ti,o would be close to the leading skeleton scale-shift to,o- In this case the single scale 
setting procedure p, could give similar results. The same holds in the skeleton 



scheme, if ^2,2 is close to /5^2- This can be deduced from eq. ( |105| ) which gives, 

^2,2 - P2,2 = s, {ri'^ - 2ri')) = 2si (ti,o - to,o) , (109) 



where in the last step we used the leading order results for the scale-shifts in eq. 
and (^). In this respect it is interesting to note that applying multi-scale BLM in 
MS, one in general obtains large values for the tifl scale-shift since = 3.385 is not 
close to (32,2 of the physical effective charges. For example, when applying BLM to 
'^-D(aMs) one obtains fco,o = 0.707(5 and ki^ = 0.366 10~^Q. This can be contrasted, 
for instance, with the BLM scales for aoiav)'- ^0,0 = 1-628(5 and ki^o = 2.487(5- 
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